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Summary. We show that a nonlinear dynamical system in Poincare-Dulac 
normal form (in R") can be seen as a constrained linear system; the constraints 
are given by the resonance conditions satisfied by the spectrum of (the linear 
part of) the system and identify a naturally invariant manifold for the flow of the 
"parent" linear system. The parent system is finite dimensional if the spectrum 
satisfies only a finite number of resonance conditions, as implied e.g. by the 
Poincare condition. In this case our result can be used to integrate resonant 
normal forms, and sheds light on the geometry behind the classical integration 
method of Horn, Lyapounov and Dulac. 



Introduction 

Normal forms are central to our understanding of nonlinear dynamics around 
known solutions, in more ways than we could recall here; see e.g. ^]. 

In this note we will consider dynamical systems - or, equivalently, vector 
fields - in R" which are in normal form around a regular critical point (for 
normal forms in this context, see jl], |2|, ^, [l6|], or the short introductions 
given e.g. in g | g @). 

We recall that if the linear part of the system is nonresonant (this and other 
relevant definitions will be recalled in section 1), the normal form is linear; we 
will thus assume the linear part is resonant. If the system is in normal form, all 
nonlinear terms will be resonant. 

We will show that the time evolution of resonant nonlinear terms can always 
be written as a linear automorphism of their set; this will allow to map the 
nonlinear system to a "parent" system of larger dimension, but linear. 

* Supported by "Fondazione CARIPLO per la ricerca scientifica" under project "Teo- 
ria delle perturbazioni per sistemi con simmetria" ; E-mail: gaeta@berlioz.mat.unimi.it or 
gaeta@romal.infn. it 
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If the linear part satisfies some additional conditions - in particular, the 
Poincare condition, i.e. if the spectrum of the linear part lies in a Poincare 
domain - the linear system has finite dimension, and the construction considered 
here leads to an explicit and simple method of integration. 

Resonant normal forms satisfying the Poincare condition are well known to 
be integrable using a "triangular" property of equations in resonant normal form 
which was already remarked by Dulac ]5|], and which guarantees integrability 
also with a classical procedure (Dulac attributes this to Horn and Lyapunov); 
for related modern ideas see Q. The present note shows how this very clas- 
sical result is related to the possibility of writing a resonant normal form as a 
constrained linear system. 

Actually, this note implements for resonant normal forms ideas which were 
put forward (in the hamiltonian context) by several authors, in particular Kazh- 
dan, Kostant and Sternberg (ll) (see also (l^, [ll| ) . 

Their point is, roughly speaking, that while in many cases one uses a sym- 
metry quotient (Marsden-Weinstein reduction) to reduce a nonlinear system to 
a lower dimensional one, and maybe integrate it, in other cases the nonlinear 
system can be seen as the reduction to a lower dimensional invariant manifold 
of a linear system living in a higher dimensional space, and can be integrated 
by "enlargement" rather than by reduction. 

This mechanism is at work in systems of Lax type (e.g. Calogero-type sys- 
tems), where an n-dimensional nonlinear system is embedded into a higher 
dimensional - typically n 2 dimensional, resorting to [n x n) matrices - linear 
system. It was even conjectured Q that all integrable systems actually origi- 
nate in higher dimensional linear systems, which we observe only through their 
projection to a lower dimensional nonlinear manifold. 

The present work shows that the Kazdan-Kostant-Sternberg mechanism can 
work also for non-hamiltonian systems, and that it applies in particular to res- 
onant normal forms. 

After a preliminary version of this note was circulated, prof. Walcher pointed 
out that the basic ideas behind the treatment of resonant normal forms given 
here can also be found in ||l6f ; however I believe that the implementation given 
in this note is somewhat simpler, and clarifies the geometry involved in the 
problem and in the integration method. 
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1 Normal forms 



In this section we fix notation and collect several definitions, and properties of 
normal forms, to be used in the following. 

Let us consider a C°° function / : R" — > R" such that /(0) = 0, expanded 
in a series of homogeneous terms fk(x), where fk(ax) = a k+1 fk(x) for all real 
number a. This defines a dynamical system x = f(x) in Vq = R" having a criti- 
cal point at the origin; we also rewrite this, using the expansion in homogeneous 
terms and singling out the linear part, as 

oo 

x = Ax + F(x) = Ax + 53/ fc (a:) ■ C 1 ) 

fe=i 

We will assume that the matrix A := (Df)(Q) is in Jordan form. 

Remark 1. Transforming a general matrix to Jordan form can be, in prac- 
tice, very hard for large n\ see |l4|] for a discussion of this point in the context 
of normal forms theory. 

We consider several vector fields associated to the system and to the decom- 
position given in (1), to be used in the following; we always write di for djdx 1 . 
As well known the matrix A can be decomposed into a semisimple and a nilpo- 
tent part, mutually commuting, which we denote as A s and A n ; as [A s , A n ] = 0, 
we also have [A s , A] = = [A n , A}. 

We will use X f = f i {x)d i , X A = {Ax) l d u X = {A s x) l d t , and X F = 
F l (x)di; we also consider Xg = (A + x) l di; as long as we deal with real matrices 
A, this reduces to X t = (A T x) t d l . Note that [X ,X A ] = = [X ,Xt\. 

Finally, we recall that a matrix A is said to be normal if it commutes with 
its adjoint, [A, A + ] = 0; obviously this is equivalent to the condition that 
=0. In general this is not satisfied, and Jacobi identity only guar- 
antees that [X , [Xa, X(}] = 0. 

Let (Ai, A„) be the eigenvalues of A (we denote by a their ensemble, i.e. 
the spectrum of A), and let (ei, e„) be the basis in Vo = R" with respect to 
which the x % coordinates are defined. We will use the multiindex notation 

t m ^Ml ,Ce 

. \ h n 

Then the vector v^) :— x^e a is resonant with A (or resonant for short) if 

n n 

(p - A) := /J,i Xi = A Q with pn > , \fx\ := ^A*i > 2 . (2) 

1=1 2 = 1 

Note that when we define and determine resonant monomials and vectors, we 
can as well consider A s rather than A. 

The space of vectors resonant with (the semisimple part of) A is defined as 
the linear span of such vectors; we will consider a basis (vi, v r ) in this. Thus 
F is resonant if and only if F = ciWi for some constants Cj. 
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Remark 2. We stress that r could be infinite, but will always be finite if 
A admits a finite number of resonance relations, e.g. if a satisfies the Poincare 
condition. On the other side, if the \ satisfy a "master resonance" relation, i.e. 
Y]. fiiXi = with fii as in (2) and > (notice in this case a cannot satisfy 
the Poincare condition), then there is an infinite number of resonances. This is 
in particular the case if there is a zero eigenvalue or if we deal with a resonant 
Hamiltonian system. 

A monomial x^ such that (2) is satisfied for some r £ (1, n) is called a res- 
onant monomial. We consider a basis of resonant monomials {^(x), ...,4> r (x)}, 
and their linear span; this is a linear vector space V\ (in the space of scalar poly- 
nomials on R"). We choose a basis {pi, ...,p r } in this. Here p l corresponds to 
4> l (x), i.e. the scalar polynomial 53i=i c i4> l ( x ) is represented in V\ by the vector 

We say (see e.g. ||) that (f ) is in Poincare-Dulac normal form if the vector 
fields Xi and Xp commute: 



This implies that all nonlinear terms are resonant with A (i.e. with A s ); however 
not all resonant terms will satisfy (3) when A n ^ 0, see e.g. example 3 below. 

Notice that in general [X^A^] ^ 0: thus, unless A is normal, we cannot 
affirm that Xg (or Xa) commutes with Xf. 

However it is easy to see from (3) that for systems (1) in Poincare-Dulac 
normal form, both Xa and Xp commute with Xq, and therefore 



When the system satisfies (4) - although (3) is possibly not satisfied - i.e. if 
F is resonant with A s , we say that it is in normal form with respect to ^ or, 
that it is in seminormal form. 

As well known, starting from any dynamical system (or vector field) of the 
form (1), we can arrive to a dynamical system (or vector field) in Poincare- 
Dulac normal form by means of a sequence (in general, infinite) of near-identity 
transformation obtained by means of the Poincare algorithm; these combine 
into a near-identity transformation H defined by a series which is in general 
only formal. The same applies for seminormal forms. 

However, one can guarantee the convergence of the series on the basis of 
properties of the spectrum a of the linear part A. In particular, it was already 
known to Poincare and Dulac that convergence is guaranteed if the convex hull 
of a in the complex plane does not include the origin; in this case we say that 
a belongs to a Poincare domain, or that A satisfies the Poincare condition, or 
also that a is a Poincare spectrum. This relevant property makes the study of 
normal forms with Poincare spectrum specially interesting. 

It is easy to see - and again was well known to Poincare and Dulac - that 
if er satisfies the Poincare condition, then only a finite number of resonances 




. 



(3) 



[X ,Xf] = 



(4) 
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is present, i.e. the Poincare-Dulac normal form is finite (conditions ensuring 
finiteness of the normal form are discussed in R] ) . 

Remark 3. Other conditions guaranteeing convergence of H, on the basis 
of a and of symmetry properties of the normal form, are also known, see e.g. 
the review 0] ; in such cases one is not guaranteed to have an integrable normal 
form and is thus no surprise that the method presented here does in general not 
apply. 

2 Normal forms as reduction of a linear system 

In this section we will give a very elementary procedure to associate to a non- 
linear system (1) in resonant normal form a linear system £ — Bt; in a real 
vector space V ~ V © V\ (this can be seen as a trivial bundle tt : V — ► Vo 
over Vo) and a set of constraints £ so that £ = identifies a smooth (actually, 
algebraic) submanifold M. C V, invariant under the flow of the linear system. 
The resonant normal form is just the reduction of the linear system to M.. 

If the normal form is finite then V — R w for some finite N > n, and 
this procedure also provides a way to explicitely and elementarily integrate the 
system in normal form, as discussed in the section 3. 

2.1 General construction 

Let Vi, 4>i(x) and pi be as defined above; we assume for ease of language that 
r is finite (as mentioned above, this is the case for a satisfying the Poincare 
condition) . We consider the real vector space V = Vo ® Vi = R n+r with basis 
vectors (ei, e„; pi, p r ) and coordinates £ = (x 1 , ...,x n ;w 1 , ...,w r ). 

We consider a "parent" dynamical system £ = [a vector field — 

^(Oi^/dS,^) ] in V defined as follows: first we rewrite (1) substituting w 1 , ...,w r 
for ^(x), <jf{x) (this gives the evolution equation for the x's); then we assign 
time evolution for the w by dw l /dt = (d<j) l (x)/dx^)(dx : ' /dt). Having written 
these equations, we will now consider the x and w as indipcndcnt quantities. 

It is clear that, by construction, the manifold Ai C V defined by 

£* := w i -4> i (x 1 ,...,x n )=0 Vi = l,...,r (5) 

is invariant under the flow of the system we have defined in this way, i.e. : 
M -» TM. 

The M. defined by (5) is an algebraic manifold (since the <f> are polynomials) 
and is tangent in the origin to the linear space R™ defined by w l — 0, i.e. to 
Vq C V (since the (f> are nonlinear functions of the x). It is also obvious from 
(5) that it is a global section of tt : V — > Vq. Thus a smooth dynamics on A4 
projects globally to a smooth dynamics in Vq. 

On the invariant manifold A4, the parent system is equivalent by construc- 
tion to the original one. 
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Remark 4. In the fiber bundle notation, and denoting by ir* the lift of 
7T : V ->• V to 7T* : TV -> TVb and by Y^ 1 the restriction of F v , to X, we have 
7 r*[Y 1 ^]=X / . 

We will show in the next subsection that the evolution equation we obtain 
for £ := (x, w) is linear, i.e. £ = W(^) reduces to £ = 2?£, with B a matrix. 

The evolution equations constructed in this way have also several other gen- 
eral features in common; we mention these omitting the elementary proof. 

(1) We actually obtain a "block triangular" evolution equation: the evolution 
of the w depends only on the w themselves, while that of the x depends on the 
x and the w together. 

(2) The system of ODEs we obtain is also, if the coordinates are properly 
arranged]^ triangular in proper sense. 

(3) The eigenvalues of B (see point 1) are given by (Ai, A„; A Q1 , X ar ), 
where the A,'s are the eigenvalues of A, and on is the a associated to the resonant 
monomials <pi, see (2); thus, we always have multiple eigenvalues. 

2.2 Proof of the linearity of the parent system 

We will now substantiate the assertion that the parent system obtained accord- 
ing to the procedure given above will always be linear. We will actually proof 
this in two ways, i.e. both algebraically and geometrically. 

Algebraic proof. 

Summation over repeated indices will be understood, and we will denote by 
v(i) the multiindex such that Vj = Stj. 

Let us consider a resonant monomial w = , and say it satisfies (A • /i) := 
Xi/ii — X a . We denote by W a ;i;a-, i = l,--;<l(ct) the resonant monomials x a 
such that (A • a) = X a (with the same fixed a). Then we have that (for / in 
seminormal form) under x = f{x) the time evolution of w is given by Xf(w), 
i.e. 

where we have used the decomposition (1) and the fact all the nonlinear terms 
must be resonant. 

We assume that A is in Jordan form, so that A s — diag(Ai, A n ), and 
{A n ) % - — rjij is different from zero (and equal to one) if and only if j = i + 1 and 
x l , belong to the same Jordan block; this implies of course that A, = Xj. 

Notice that terms with \ii = are absent from the sum over i; we can 
therefore assume fXi ^ 0. Under this condition, and using the assumption that 
A is in Jordan form, we can rewrite 

w = (A • n) x? + m rjij X »-»W+»U) + Ca ^ x ^-^)+<y . ( 6 ) 

1 in general, we can have to pass to complex coordinates in the x space and to have a change 
of variables in the w space to obtain this, see example 4 
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The first term on the r.h.s. is nothing else that X a w. We want to check that 
the other terms are also (the sum of) resonant monomials; in order to do this 
we do not have to worry about the scalar coefficients in fronts of them. 

The monomials appearing in the second term of the r.h.s. are of the form 
x'-f = a^-^OHKfc) j an d wc can assume \Xi ^ and rjij ^ (or the corresponding 
monomial would not be present in the sum). For these we have 

(A • ip) := \ l Lp l = (A • fj,) - Aj + Xj = (A • fi) = X a ; 

we have used the fact that rj^ ^ implies Aj = Xj, and the resonance relation 
satisfied by w itself. Thus the second term in the r.h.s. of (5) is the sum of 
resonant monomials (with the same aasai). 

The monomials appearing in the third term are of the form x v = x IJ, ~ v ^ >+ ' 7 , 
where (A • a) = Aj and we can assume ^ ^ (or the corresponding monomial 
would be absent from the sum). We have now 

(A • tp) := Xufi = (A • fi) - Xi + (A • a) = X a ; 

again we have a sum of resonant monomials (with the same a as to). 

This concludes the proof that the right hand side of (6) can be written as 
a linear combination of resonant monomials, i.e. that the evolution equation 
constructed according to our procedure is linear. A 

Remark 5. Notice that we have actually proved something more, i.e. that 
if w = x^ with (A • fj,) = X a , only resonant monomials w = x v with (A • ir) = X a 
(with the same a as above) will appear in this linear combination. That is, 
the matrix B will be a block one, where the blocks correspond to resonant 
monomials identified as described here. 

Geometric proof 

A more geometric (but equivalent) proof could be obtained by considering a 
basis of (nonlinear) resonant vectors x^ek (k = 1, ...,n) and the corresponding 
vector fields = {x^)dk {i = 1, see above). These obviously generate a 
Lie algebra Q (recall that the resonance condition is equivalent to commutation 
with Xq, and notice that the commutator of two vectors built from nonlinear 
monomial terms will never be a linear vector field). 

One should then check that Q is globally invariant under time evolution, i.e. 

that [x f ,g\ C g. 

Take £ Q: we have [Xf,X$] — [Xa, X^] + [Xp, X^], and the second term 
is by definition in Q. To see that the first is also in Q, we have to check the 
vanishing of [Xq, [X a, X ,/>]]; using the Jacobi identity, this is [Xa, [Xo,X$]] — 
[X^, [X ,X A ]], and both terms vanish separately. The proof is complete. A 

Remark 6. Notice that if A is not normal, the set of resonant monomials 
defining vectors in normal form with respect to A + would not, in general, be 
closed under time evolution. This is also immediately seen from the alternative 
geometrical proof by remarking that if we substitute Xi for Xq, we have (for A 
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not normal) [X^X^] ^ 0. Thus for non- normal A we cannot limit to consider 
vectors in normal form, but have to consider the set of all resonant vectors. This 
will also be clearly shown in example 3 below. 



2.3 Truncated normal forms 

When the normalizing transformation is not convergent, the normal form is not 
conjugated to the original system. However, in such a case one can consider 
normalization only up to a sufficiently low degree N (in practice this is deter- 
mined by either the computational limits or the optimal degree on the basis 
of convergence properties of the truncated series); in this way one obtains a 
system which is in normal form up to order N . This system can be truncated 
at order N - thus obtaining a truncated normal form - and the relation of such 
a truncation with the full system will then be studied via other techniques in 
perturbation theory Jl], ||, [|, 
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The block structure of the B matrix, as determined above (see in particular 
remark 5), explains when this truncation will result in a closed parent system. 

Indeed, consider all the resonance relations (2); let m_(a) and m+(a) be 
the smaller and greater values of |/i for which a relation (fx- X) = X a is satisfied. 
It follows from remark 5 that if 



N g [m-(a),m+(a)] Va = l,...,n , (7) 

then the truncated normal form at order N is mapped, by the procedure dis- 
cussed in this note, to a closed linear system (note that for a finite dimensional 
system we always have a finite number of resonances with \fi\ < N, for any finite 

AO- 

Thus the integration procedure discussed here is also of use in cases where 
the normalizing transformation is not convergent in any neighbourhood of the 
origin. 

On the other hand, it should be stressed that if the system admits an infinite 
number of resonances, then (7) can be satisfied only for N < m_(a) for all a, 
but in such case the truncated normal form is trivial (it reduces to the linear 
part of the system). Thus, the truncated normal form will not result in a closed 
finite dimensional linear parent system, see example 4 below. Such a situation 
is met when a "master resonance" is present (see remark 2), and in particular 
when dealing with resonant hamiltonian systems. 



3 Integration of normal forms 

The strategy to integrate normal forms via the parent linear system £ = _B£ 
we have defined is rather obvious: this rests on the dynamical invariance of the 
manifold defined by (5) and can be divided into three steps. That is, 

1. Step 1. For £ = (x; w), determine the general solution of the linear equa- 
tion £ = B£ in V = R™ +r , say with solution £(£) where £(0) = £o = 
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(xo,Wq) is the initial datum. This will depend on the n + r arbitrary 
constants (xo,wq). 

2. Step 2. Restrict the general solution to the invariant n-dimensional sub- 
manifold M. C R n+r defined by w i — ipiix 1 , ...,x n ). This will depend on 
the n arbitrary constants Xq . 

3. Step 3. Project the general solution (x(t), w(t)) on Ai C V to the subspace 
Vq = R™ spanned by the x variables, i.e. extract x(t) forgetting about 
w(t). 

As discussed above, the correspondence between the original nonlinear sys- 
tem and the restriction of the parent system to the invariant manifold M is 
guaranteed by construction, and projection is globally well defined as M. is 
identified by the algebraic equations (5). It is therefore clear that this pro- 
cedure will indeed provide the most general solution to the original nonlinear 
system in Vq. 

This strategy will be particularly simple, and successful, when there is only 
a finite number of resonances, and in particular when a belongs to a Poincare 
domain. 

It should be stressed that if a belongs to a Poincare domain, the normal 
forms could of course also be integrated directly: indeed the corresponding 
system is nonlinear but, as remarked by Dulac Q, always in triangular form 
(cf. the properties of the parent system mentioned at the end of subsection 
2.1). Namely, we can always write x 1 — A^x^ + <I> 4 (x) in such a way that 
d^/dx^ — for j > i. It is then possible to solve the equations recursively, 
starting from the linear one for x 1 ^) and having at each step a linear equation 
with a forcing term which is a nonlinear but explicitely known function of t. 

The procedure proposed here is equivalent to the one considered by Dulac 
(and attributed by him to Horn and Lyapunov) from the analytic point of viewQ, 
but it has the advantage of showing how the nonlinear (normal form) system is 
obtained by restrictrion (on the submanifold M.) of a linear system via nonlinear 
constraints, clarifying the geometry involved in the integrability and integration 
of resonant normal forms with Poincare spectrum, and the connection with 
topics in modern integrable systems theory (ll], It also shows that these 
ideas, developed in the hamiltonian context, can be fruitfully applied to more 
general dynamical systems. 

4 Examples 

Systems in normal form are specially interesting if a small neighbourhood of 
the origin is dynamically invariant, i.e. if the critical point is stable (in this 

2 Note it can be more convenient computationally, as it only requires to solve linear systems; 
in particular it will be conveniently implemented on computers via algebraic manipulation 
languages. 
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case the evolution will remain in the domain of analyticity of the normalizing 
transformation); thus we are mainly interested in cases where the real part of the 
eigenvalues is negative (or zero). However, for ease of notation we will consider 
examples with positive eigenvalues; the stable situation is recovered by a time 
reversal. Also for ease of notation, we will write all vector indices as lower ones; 
the Cj will be arbitrary real constants. 



Example 1. 

Let us consider n = 2, with coordinates (x, y) and 



A = 



1 

k 



with k a positive integer; notice here a = {1, k} is in a Poincare domain. There 
is only one resonance fj, — (k, 0) (with a ~ 2), and the only resonant monomial 
is (j>(x, y) = x k . The Poincare-Dulac normal form is 

x = x 

y = ky + cix k 

with ci an arbitrary coefficient. 

Thus, following our procedure, we set w = k(x k ~ 1 )x = kw; the system 
obtained in V is 

x = x 

y = ky + ciw 

w = kw 



and the constraint £ is given by w — x k =0. 
The solution to the system is given by 

x(t) = x e* , y(t) = y e kt + ( Cl kw )te kt , w(t) = w e kt ; 

obviously the submanifold M identified by w — x k is invariant under this flow, 
and the projection of solutions on M to R 2 = (a;, y) is simply 

x(t) = x e* , y(t) = [y + { Cl kx k a )t]e kt 



Example 2. 

Let us consider n = 3, with coordinates (x, y, z) and 

/l 0\ 
A = 2 ; 
\0 5/ 

notice here a — {1, 2, 5} is in a Poincare domain. There are four resonances: 





(2,0,0) 


(with a 


— 2 and \fj,\ 


= 2) 




(1,2,0) 


(with a 


— 3 and \/j,\ 


= 3) 


M = 


(3,1,0) 


(with a 


— 3 and 


= 4) 


1 n = 


(5,0,0) 


(with a 


— 3 and \n\ 


= 5) 
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and correspondingly we have 



(pi= x 2 , (f> 2 = xy 2 , (f> 3 = x 3 y , 04 = x 5 . 



The normal form is written as 



x = x 



y = 2y + cix 2 

z = 5z + c 2 xy 2 + c 3 x 3 y + c 4 x 5 



with Cj arbitrary real coefficients. 

By our procedure, we obtain a system in V = R 7 , given by 



X = 


X 


y = 


2y + awi 


z = 


5z + c 2 w 2 + c 3 w 3 + C4W4 


W\ = 




U>2 = 


5w 2 + 2ciw 3 


W 3 = 


5w 3 + ciw 5 


W4 = 


5u>4 



(notice this has the block structure discussed in section 2). The general solution 
to this (writing w,(0) = 7 j for ease of notation) is 

x(t) = x e* , y(t) = (y a + jit)e 2t 

z(t) = [z + (72C2 + 73C3 + 74C 4 )i + (l/2)(2 7 3Cic 2 + 7 4 c 3 )t 2 ] e 5t 
w 1 ^) = 7 ie 2t , w 2 {t) = ( 72 + 2c ll3 t)e 5t 
w 3 {t) = ( 73 + 74 i)e 5t , w 4 (t) = 74 e 5t 

which once restricted to the manifold M (here identified by £\ := w\ — x 2 , 
£ 2 := w 2 — xy 2 , £ 3 := w 3 — x 3 y, £4 := w 4 — x 5 ) and projected to R™, gives 

x(t) = x e* 

y(t) = (y + x 2 t) e 2t 

z(t) = [z + (c 2 x y 2 + c 3 xlyv + c 4 a^) t + (20^x^0 + c 3 x%) (t 2 /2)] e 5t 

Example 3. 

Let us consider n = 3, with coordinates (x, y, z) and 



notice here a — {1, 1, 2} is in a Poincare domain, and A n ^ for 77 ^ 0. 
There are three resonances, all of them with a — 3 and \fi\ = 2, i.e. 




M= (2,0,0) , fi= (1,1,0) , n= (0,2,0) 
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and correspondingly 

01 = x 2 , 2 = xy , 03 = y 2 . 
The seminormal form with respect to the semisimple part A s of A is thus 

x = x + rjx 

y= y 

z = 2z + C\X + c 2 xy + c 3 

The normal form with respect to the full A is the same for 77 = (which means 
A = A s ), and for 77 ^ is obtained setting C2 = c 3 = 0, as readily seen by 
considering Xp — (c\x 2 + c 2 xy + c$y )d z and imposing [Xt,Xp] = 0. 
The associated linear system in V = R 6 is £ = B£, with 



B = 



( l 
















V 


1 




















2 


Cl 


C2 


C3 











2 


V 

















2 


77 


Vo 














2/ 



for the seminormal form, and the same - with C2 = C3 = if 77 is nonzero - 
when we normalize with respect to the full A. 

Writing again Wi(0) = 7^, the general solution to this linear system is 

x(t) = Xoe* 

y{t) = (2/0 + VXotfe* 

z(t) = [zq + (C171 + c 2 72 + c 3 j 3 )t + (C172 + r]c 2 j3)(t 2 /2) + 77Ci7 3 (t 3 /6)] e 2t 
wi(t) = (7i + 72* + ?773t 2 /2)e 2t 
w 2 (t) = (l2+Vl3t)e 2t 
w 3 (t) = l3,e 2t 

Restricting to M. (which in this case is identified by £\ := w\ — x 2 , £2 := W2 — xy, 
£3 '■= W3 — y 2 ) and projecting to R 3 gives 

x(t) = xq e* 

y(t) = (2/0 + V x o t) e* 

z(t) = [z + (ciXq + c 2 x y + c 3 y 2 ) t + 

+ { Cl x y + vc 2 y 2 ) (t 2 /2) + V c iy 2 (i 3 /6)] e 2t 

If 77 7^ and the system is in normal form with respect to the full matrix A 
(so that C2 = c 3 = 0) then z(t) simplifies to 

z(t) = [zo + c 1 x 2 t + c 1 xoyo(t 2 /2) + r 1 c 1 y 2 (t 3 /6)]e 2t . 

Notice that (for 77 ^ 0) the time evolution of 0i depends on <f>2, and through 
this on fa as well; thus, it would not be possible to consider a parent system 
involving only terms in normal form with respect to A + . For r\ — the equations 
for the Wi decouple, but then all the 0j are allowed in the normal form. 
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Example 4. 

Let us consider n = 2 with coordinates (x, y) and 




Here a = does not satisfy the Poincare condition. We have a master 

resonance Ai + A 2 = 0, and hence an infinite number of resonances, given by 

f p = (k + 1, k) i.e. (k + l)Ai + £;A 2 = Ai , a = 1 
\ A* = (*:,* + 1) i-e- fcAi + (k + 1)A 2 = A 2 , a = 2 

and the resonant monomials are given by 

fam-i = (x 2 + y 2 ) m x = p m x , 2m = (x 2 + y 2 ) m y = p m y , 

where p — (x 2 + y 2 ). For ease of notation we introduce coordinates q m = 

W2m-1 — 02m- 1) Pm = W 2m = (p2m- 

The normal form is then 

y= +x + EZA x2 + y 2 )HbkX + a k y) 
and this is rewritten following our procedure as an infinite linear system: 

x= -y + YX=ii a k1k-bkPk) 

y= + X + Y.k=i ( b k<lk + a k Pk) 

<lk = -Pk + 12T=i i( 2k + l)a s qk+ s - b s pk+s] 

Pk = +Qk + J27=i [bsQk+s + (2fc + l)a s pk +s ] 

Notice that truncating this at order N — 2k+l, i.e. projecting it down to the 
linear subspace of (a;, y) and the (p m , q m ) with m < k, we get a finite dimensional 
linear system (all sums go then from 1 to k). However, the dynamics of this 
projected system can and in general does fail to reproduce even qualitatively 
the full dynamics; actually, it will also not preserve the manifold M. To see 
this, it suffices to consider the truncation at k = 1, N = 3, which reads 

x= -y + (aigi - hpi) 
y= +x + (hqx + a\p{] 
Qi = -Pi 

Pi = +qi ■ 

Example 5. 

Let us briefly consider the situation mentioned in subsection 2.3; let n = 4 with 
coordinates Xi (i = 1, 4) and 

A = diag (1,2,3,10) ; 
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here again a is a Poincare spectrum. We want to consider normal forms trun- 
; the resonances of order < 3 are given by 

V =(2,0,0,0) (with a = 2 and |/z| = 2) 
iu= (1,1,0,0) (with a = 3 and j/zj = 2) 
/z= (3,0,0,0) (with a = 3 and j/zj = 3) 

There are then several resonances of higher order, with 4 < |/z| < 10, all of them 
with a = 4. 

We will introduce 

4>i = x\ , 02 = 2:12:2 , <fe = ar? . 
The normal form corresponding to this linear part A is 

X\ = X\ 

x 2 = 2x 2 + c\x\ 

is, = 3x 3 + C2X1X2 + C3X1 

±1 = 10x 4 + 0(|x| 4 ) 

where the Ci are real constants, and the truncated normal form of order three 
is obtained by dropping the term 0(|x| 4 ) in the above. 

By applying our procedure to this truncated normal form we get 

X\ — X\ 

x 2 = 2x 2 + Ci4>t 

£3 = 3x 3 + C 2 (p2 + C 3 3 

X\ = 10x 4 

(pi = 2 (pi 

<p2 — 3 (p2 + Ci 03 
<Pl — 3 03 

We stress that this is not a truncation of the obtained equations: once we deal 
with the truncated normal form, we obtain exactly this closed form (indeed, no 
(p depends on 14), as shown in subsection 2.3. 
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